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Computing fundamental domains 
for Fuchsian groups 

par John Voight 



Resume. Nous exposons un algorithme pour calculer un domainc 
de Dirichlet pour un Fuchsian groupe T avec aire cofinis. En 
consequence, nous calculons les invariants de T et une presentation 
cxplicite finis pour T. 

Abstract. We exhibit an algorithm to compute a Dirichlet do- 
main for a Fuchsian group T with cofinite area. As a consequence, 
we compute the invariants of T, including an explicit finite pre- 
sentation for r. 

Let r C PSL2(1R) be a Fuchsian group, a discrete group of orientation- 
preserving isometries of the upper half-plane fj with hyperbolic metric d. 
A fundamental domain for T is a closed domain D C f) such that: 

(i) YD = and 

(ii) gD° n D° = for all g G T \ {1}, where ° denotes the interior. 
Assume further that T has cofinite area, i.e., the coset space X = has 
finite hyperbolic area n{X) < oo; then it follows that V is finitely generated. 

In this article, we exhibit an algorithm to compute a fundamental domain 
for V; we assume that T is specified by a finite set of generators G C SL^K) 
with K 1R n Q a number field, and we call V exact. Suppose that p G fj 
has trivial stabilizer T p = {1}. Then the set 

D{p) = {z G fj : d(z,p) < d(gz,p) for all g G V}, 

known as a Dirichlet domain, is a hyperbolically convex fundamental do- 
main for r. The boundary of D(p) consists of finitely many geodesic seg- 
ments or sides. We specify D(p) by a sequence of vertices, oriented coun- 
terclockwise around p. The domain D(p) has a natural side pairing: For 
each side s of D(p), there exists a unique side s* and g G T \ {1} such that 
s* = gs, and the set of such g comprises a set of generators for T. 
Our main theorem is as follows. 

Theorem. There exists an algorithm which, given an exact Fuchsian group 
T with cofinite area and a point p G with T p = {1}, returns the Dirichlet 
domain D{p), a side pairing for D{p), and a finite presentation for V with 
a minimal set of generators. 
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This algorithm also provides a solution to the word problem for the 
computed presentation of T. 

Of particular and relevant interest is the class of arithmetic Fuchsian 
groups, those groups commensurable with the group of units 0\ of reduced 
norm 1 in a maximal order O of a quaternion algebra B defined over a 
totally real field and split at exactly one real place. Alsina-Bayer [1] and 
Kohel-Verrill [18] give several examples of fundamental domains for arith- 
metic Fuchsian groups with F = Q. Our work generalizes that of Johansson 
[15] . who first made use of a Dirichlet domain for algorithmic purposes: he 
restricts to the case of arithmetic Fuchsian groups, and we improve on his 
methods in several respects (see the discussion preceding Algorithm 1 2 . 5 1 and 
the reduction algorithms in §4). 

The algorithm described in the above theorem has the following appli- 
cations. The first is a noncommutative generalization of the problem of 
computing generators for the unit group of a number field. 

Corollary. There exists an algorithm which, given an order O C B of a 
quaternion algebra B defined over a totally real field and split at exactly 
one real place, returns a finite presentation for 0\ with a minimal set of 
generators. 

We may also use the presentation for T to compute invariants. The group 
r has finitely many orbits with nontrivial stabilizer, known as elliptic cycles 
or parabolic cycles according as the stabilizer is finite or infinite. The coset 
space X = r \ fj can be given the structure of a Riemann surface, and 
we say that T has signature (g; mi, . . . ,mt,s) if X has genus g and Y has 
exactly t elliptic cycles of orders mi, . . . ,mt £ Z>2 and s parabolic cycles. 

Corollary. There exists an algorithm which, given T, returns the signature 
ofT and a set of representatives for the elliptic and parabolic cycles in T. 

Finally, we mention a corollary which is useful for the evaluation of 
automorphic forms. 

Corollary. There exists an algorithm which, given T and z,p £ 5) with 
T p = {1}, returns a point z' £ D{p) and g £T such that z' = g(z). 

The article is organized as follows. We begin by fixing notation and 
discussing the necessary background from the theory of Fuchsian groups 
(§1-2). We then treat arithmetic Fuchsian groups and give methods for 
enumerating "small" elements of the group 0\, with OcBa quaternion 
order as above (§3). Next, we describe the basic algorithm to reduce an 
element g £ T with respect to a finite set G C T (§4). We then prove the 
main theorem (§5) and conclude by giving two examples (§6). 

The author would like to thank the Magma group at the University 
of Sydney for their hospitality, Steve Donnelly and David Kohel for their 
helpful input, and Stefan Lemurell for his careful reading of the paper. 
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1. Fuchsian groups 

In this section, we present the relevant background from the theory of 
Fuchsian groups; suggested references include Katok [HI Chapters 3-4] 
and Beardon [2^ Chapter 9]. Throughout, we let T C PSL^R) denote a 
Fuchsian group with cofinite area, which is finitely generated by a result of 
Siegel [161 Theorem 4.1.1], |12[ §1]. To simplify, we will identify a matrix 
g G SL 2 (M) with its image in PSL 2 (M). 

Throughout this section, let p G be a point with trivial stabilizer 
T p = {1}. Almost all points p satisfy this property: there exist only finitely 
many p with T p ^ {1} in any compact subdomain of fj, and in particular, 
the set of p G fj with T p ^ {1} have area zero. In practice, with probability 

1 a "random" choice of p will suffice. 

We define the Dirichlet domain centered at p to be 

D(p) = {z G fj : d(z,p) < d(gz,p) for all g G T}. 

The set D(p) is a fundamental domain for T, and is a hyperbolic polygon. 
More generally, we define a generalized hyperbolic polygon to be a closed, 
connected, and hyperbolically convex domain whose boundary consists of 
finitely many geodesic segments, called sides, so that a hyperbolic polygon 
is a generalized hyperbolic polygon with finite area. 

Let D C fj be a hyperbolic polygon. Let S = S(D) denote the set of 
sides of D, with the following convention: if g G T is an element of order 

2 which fixes a side s of D, and s contains the fixed point of g, we instead 
consider s to be the union of two sides meeting at the fixed point of g. We 
define a labeled equivalence relation on S by 

P = {(g,s,s*) : s* = g(s)} CT x (S x S). 

We say that P is a side pairing for D if P induces a partition of S into 
pairs, and we denote by G{P) the projection of P to V. 

Proposition 1.1. The Dirichlet domain D(p) has a side pairing P, and 
the set G(P) generates T. Conversely, let D C be a hyperbolic polygon, 
and let P be a side pairing for D. Then D is a fundamental domain for 
the group generated by G(P). 

Proof. The first statement is well-known [21 Theorem 9.3.3], |16^ Theorem 
3.5.4]. For the second statement, we refer to Beardon [21 Theorem 9.8.4] 
and the accompanying exercises: the condition that /u(Z?) < oo ensures 
that any vertex which lies on the circle at infinity is fixed by a hyperbolic 
element [121 §1]- □ 

Remark 1.2. The second statement of Proposition 11.11 extends to a larger 
class of polygons (see J2J §9.8]), and therefore conceivably our results extend 
to the class of finitely generated non-elementary Fuchsian groups of the first 
kind. For simplicity, we restrict to the case of groups with cofinite area. 
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We can define an analogous equivalence relation on the set of vertices of 
D, and we say that a vertex v of D is paired if each side s containing v is 
paired to a side s* via an element g G G such that gv is a vertex of D. 

We now consider the corresponding notions in the hyperbolic unit disc 
2), which will prove more convenient for algorithmic purposes. The maps 

(1.1) z-p pw-p 

z — p w — 1 

define a conformal equivalence between fj and 2) with p cf>(p) = 0. Via 
the map <f>, the group T acts on 2) as 

= ^(j)- 1 C PSU(1, 1) = j± ^ ^ € PSL 2 (C) : a = d, b = c j . 

We may analogously define a Dirichlet domain D(q) for g 6 3 with T q = 
{1}, and we have 4>(D(p)) = -D(O) C 2). To ease notation, we identify V 
with by g i— > = <\>g4>~^ when no confusion can result. 

Any matrix = ^ ^\ G SU(1, 1) acts on 2), multiplying lengths 

by |</(2)| = |cz + c?|~ 2 , and therefore Euclidean lengths (and areas) are 
preserved if and only if \cz + d\ = 1. We define the isometric circle of 5 to 
be 

1(g) ={zeC:\cz + d\= 1}; 

if c / 0, then 1(g) is a circle with radius l/\c\ and center —d/c, and if c = 
then 1(5) = C. We denote by 

int(/( 5 )) = {z G C : \cz + d\ < 1}, ext(/(c/)) = {z G C : \cz + d\ > 1} 

the interior and exterior of 1(g), respectively. 

With these notations, we now find the following alternative description 
of the Dirichlet domain D(0) C 2). 

Proposition 1.3. 

(a) The domain D(0) is the closure in 2) of 

n ext (%))- 

<yer\{i} 



(b) For any g 6 SU(1, 1), we /iai>e 



d(z,0) < 



< 



> 



d(gz, 0) according as < 



r zeext(I(g)), 

z G /(5), 
^ G int(/( 5 )). 
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c d 
a 2 + b 2 + c 2 + d 2 



Proof. See Katok [16^ Theorem 3.3.5]; we note that if g G T and c = 0, then 
q = is a fixed point of g, so by hypothesis g = 1, and hence ext(/(g)) ^ 
for all g ^ 1. In particular, since T has cofinite area we note that the 
intersection in (a) is nonempty. □ 

Corollary 1.4. For any g G SU(1, 1), we have gl(g) = I(g~ l ). 

Proof. By Proposition 11.3( b). we have 

w = gz G I{g~ l ) & d(g~ 1 w, 0) = d(w, 0) 44> cf(z, 0) = d(gz, 0) 4^ z £ 1(g) 
and the result follows. □ 

Remark 1.5. One can similarly define isometric circles 1(g) for g G PSL2(M) 
acting on Sj. One warning is due, however: although _1 (D(O)) = D(p) C 
fj is again a Dirichlet domain, its sides need not be contained in isometric 
circles (as the map eft is a hyperbolic isometry, whereas isometric circles are 
defined by a Euclidean condition). Instead, we see easily that 

cTW) = U G Si : d(z,p) = d(gz,p)}, 

i.e., the isometric circle I(g^) corresponds in fj to the perpendicular bi- 
sector of the geodesic between p and g(p). In particular, if p = i then 

a somewhat lengthy calculation reveals that for g = ( ^ ^ G SL2 

this perpendicular bisector is the half-circle of radius , , , 

(a z + <r — \y 

ab + cd 

centered at — G R. 

a z + <r — 1 

The domain D(0) is also known as a Ford domain, since Proposition 
11.31 is originally attributed to Ford |1 Ij. Theorem 7, §20]. The heart of 
our algorithm (as provided in the main theorem) will be to algorithmically 
construct a Ford domain. 

2. Algorithms for the upper half-plane and unit disc 

We represent points p G Sj,D using exact complex arithmetic: see Pour- 
El-Richards [19], Weihrauch [24] for theoretical foundations (the subject 
of computable analysis) and e.g. Boehm [3], Gowland-Lester [13] for a dis- 
cussion of practical implementations. Alternatively, our algorithms can be 
interpreted using fixed and sufficiently large precision; even though one can- 
not predict in advance the precision required to guarantee correct output, 
it is likely that an error due to round-off will only very rarely occur in prac- 
tice; see also Remark 12 .61 The induced action on T) has r «-> C SU(1, 1), 
represented as matrices with exact complex entries. 

A Fuchsian group V is exact if it has a finite set of generators G C SL^-K") 
with K =— > Q n R a number field; from now on, we assume that the group V 
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is exact. Even up to conjugation in PSL2(R), not every finitely generated 
Fuchsian group is exact; our methods conceivably extend to the case where 
the set of generators G C SL2(R) are specified with (exact) real entries, 
but we will not discuss this case any further. Algorithms for efficiently 
computing with algebraic number fields are well-known (see e.g. Cohen 

We now discuss some elementary methods for working with generalized 
hyperbolic polygons in D, which are defined analogously as those in $j. 

Let 3 = {z£C: \z\ < 1} denote the closure of D and let <9£ = {z G C : 
\z\ = 1} be the circle at infinity. We represent a geodesic L in D in bits by 
four pieces of data: 

• the center c = ctr(L) 6 CU {oo}, 

• the radius r = rad(L) G R U {oo} of L, and 

• the initial point z = in(L) G D and the terminal point w G £>; 

the inital and terminal points are normalized so that the path along L 
follows a counterclockwise orientation. Although this data is redundant, 
it will be more efficient in practice to store all values rather than, say, to 
recompute c and r when needed. 

If L\,L2 C ® are geodesies which intersect at a point v G 2) \ {0}, then 
we define Z(Li,L2) to be the counterclockwise-oriented angle at v from 
the geodesies L\ to Li for the wedge directed toward the origin, so that in 
particular we have Z(L2,Z>i) = — Z(-Li,L2)- 

Example 2.1. In the following figure, we depict a geodesic and the angle 
Z(Li,L<2) 37r/8 between geodesies. 




Figure 12.11 Geodesies and angles 



We leave it to the reader to show that one can compute using elementary 
formulae the following quantities: for geodesies L±,L2, the intersection 
L\ n L2 and (if nonzero) the angle Z(Li,L2), as well as the area of a 
hyperbolic polygon. 
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Definition 2.2. Let G C T \ {1}. The exterior domain of G, denoted 
E = ext(G), is the closure in D of the set eG ext(I(g)) fl T>. 

With this definition, Proposition 11.3( a) becomes simply the statement 
that ext(r \ {1}) is the closure of D(0). 

Let G C r be a finite subset and let E = ext(G) be its exterior domain. 
Then E is a generalized hyperbolic polygon whose sides are contained in 
isometric circles 1(g) with g G G. A proper vertex of E is a point of 
intersection v G 1(g) n /(</) between two sides (with g ^ g' £ G); a vertex 
at infinity of is a point of intersection w 6 /(g) n 52) between a side and 
the circle of infinity. A vertex of E is either a proper vertex or a vertex at 
infinity. 

Definition 2.3. Let 22 = ext(G) be an exterior domain. A sequence U = 
gi , . . . , g n is a normalized boundary for E if: 

(i) E = ext(U); 

(ii) I(gi), ■ ■ ■ ,I(g n ) contain the counterclockwise consecutive sides of D; 
and 

(hi) the vertex v £ E with minimal arg(u) G (0, 27r) is either a proper 
vertex with v G I(gi) D 2(52) or a vertex at infinity with u G I(gi)- 

It is clear that for each exterior domain E, there exists a unique normal- 
ized boundary G for E: in (i) and (ii) we order exactly those gi for which 
I(gi) are sides of E and in (iii) we choose a consistent place to start. 

Example 2.4. In the following figure, we exhibit a normalized boundary 
G = {gi , <72, 53, 54}; the vertices t>i,i>2 are on the circle at infinity whereas 
^3)^4)^5 are proper. 



Figure 12.41 Normalized boundary of a generalized hyperbolic polygon 
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We now detail an algorithm which computes a normalized boundary for 
a given exterior domain. 

Algorithm 2.5. Let G C T be a finite subset. This algorithm returns the 
normalized boundary U of the exterior domain E = ext(G). 

1. Initialize 9 := 0, U := 0, and L := [0, 1]. 

2. Let 

H := {g e G : arg(I(g) H L) > 9}. 

a. If if = 0, let g G G be such that 

:= arg(in(J(s))) G + [0, 2vr) 

is minimal. 

b. If H / 0, let 5 G # be such that 

9 := arg(%)nL) G # + [0,2^) 

is minimal; if more than one such g exists, let g be the one that 
minimizes Z(L, 1(g)). 

Let U :=U U {#} and let L := n I). 

3. If U = {gi, ...,g n } and g n = g 1 , return U := {gi, . . .,g n -i}- Other- 
wise, return to Step 2. 

Proof of correctness. By definition ext(C/) is a generalized hyperbolic poly- 
gon. Suppose that E ^ ext(C7). Then there exists g £ G such that 
L = 1(g) n ext(C/) is not just a vertex of ext(J7). Consider the initial 
point z = in(L): either z lies on a side I(gi) of ext(£7) or z G 9S. 

Suppose that z G Let fj be the initial vertex of the side Sj C /(ft). 

Then in the ith iteration of Step 2 of the algorithm we have g G H, so the 
terminal vertex of Sj is proper and we are in case (b). But by assump- 
tion we have d(vi, z) < d(vi,Vi + i) since I(gi) is a geodesic, and arg increases 
along si with the distance, thus according to the stipulations of the algo- 
rithm we must have z = Uj+i. But then in order for the interior of 1(g) to 
intersect ext(C7) nontrivially, we must have Z(L,I(gi)) < Z(I(gi + i), 1(g)), 
a contradiction. 




So suppose that z G d®. Then there exists i such that z lies on the 
principal circle between the terminal point of I(gi) and the initial point of 
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I(<7j + l). But then arg(in(/(g))) < arg(in(/(<7j + i)), contradicting (a). This 
proves that (i) holds in Definition 12,31 

It is obvious that (ii) holds, and condition (iii) holds by initialization: if 
the vertex v G E with minimal arg(u) G (0, 2ir) is a vertex at infinity then 
it is found in the first iteration of the algorithm in stage (a), and if it is a 
proper vertex then it is found in the second iteration in stage (b). □ 

A Ford domain -D(O) is specified in bits by a normalized boundary G for 
D(0). We can similarly specify a Dirichlet domain D{p) by an analogously 
defined normalized boundary of perpendicular bisectors, as in Remark 1 1. 5 1 
for many purposes, it will be sufficient to represent D(p) by a sequence of 
vertices (ordered in a counterclockwise orientation around p). 

Remark 2.6. Although the intermediate computations as above are of a 
numerical sort, an algorithm to compute a Dirichlet domain accepts exact 
input and produces exact output. 

3. Element enumeration in arithmetic Fuchsian groups 

In this section, we treat arithmetic Fuchsian groups, and in particular 
we exhibit methods for enumerating "small" elements of these groups. See 
Vigneras [22] for background material and Voight [231 Chapter 4] for a 
discussion of algorithms for quaternion algebras. 

Let F be a number field with [F : Q] = n and discriminant dp. A 
quaternion algebra B over F is an F-algebra with generators a, (3 G B such 
that 



a 2 = h, P 2 = k, (3a = —aj3 

(h k\ 

with h,k G F*; such an algebra is denoted B = I — I and is specified in 



bits by h, k G F* . An element 7 G B is represented by 7 = x+ya+z(3+wa(3 
with x, y,z,w G F, and we define the reduced trace and reduced norm of 7 
by trd(7) = 2x and nrd(7) = x 2 — hy 2 — kz 2 + hkw 2 , respectively. 

Let B be a quaternion algebra over F and let TLp denote the ring of 
integers of F . An order O C B is a finitely generated Zi?-submodule with 
FO = B which is also subring; an order is maximal if it is not properly 
contained in any other order. We represent an order by a pseudobasis over 
7Lp\ see Cohen §1] for methods of computing with finitely generated 
modules over Dedekind domains using pseudobases. 

A place v of F is split or ramified according as B v = B <g>p F v = M2(F V ) 
or not, where F v denotes the completion at v. The set S of ramified places 
of B is finite and of even cardinality, and the ideal D = n^es ^|oo Pv of Zf 
is called the discriminant of B. 

Now suppose that F is a totally real field, and there is a unique split real 
place v G" S corresponding to : B M2(M). Let O C B be an order and 
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let 0\ denote the group of units of reduced norm 1 in C Then the group 
r(O) = too(^/{±l}) C PSL 2 (IR) is a Fuchsian group [HI §§5.2-5.3]. If O 
is maximal, we denote T B {1) = T(0). An arithmetic Fuchsian group Y is a 
Fuchsian group commensurable with r B (l) for some choice of B. One can, 
for instance, recover the usual modular groups in this way, taking F = Q, 
O = M 2 (Z) C M 2 (Q) = B, and T C PSL 2 (Z) a subgroup of finite index. 

An arithmetic Fuchsian group T has cofinite area; indeed, by a formula 
of Shimizu [201 Appendix], the area A = fJ,(X) = /i(r\fj) is given by 

(3-1) A = ^L-4/ 2 CF (2)d>(c,)[r B (l) : r], 

where Cf(s) denotes the Dedekind zeta function of F, and 

m = #(z F /*z F )* = N(a) [J (i - j^y) ; 

here the hyperbolic area is normalized so that 
and hence an ideal triangle has area 1/2. 

Remark 3.1. The area A is effectively computable from the formula (]3.ip . 
By the Riemann-Hurwitz formula, we have 

(3.2) A = 2g-2 + Y,e q U--\+e 00 

q ^ ' 

where e q is the number of elliptic cycles of order q £ Z> 2 in T and eoo 
the number of parabolic cycles. In particular, A £ Q; and since e g > 
implies F(£ 2(J ) i?, the denominator of A is bounded by the least common 
multiple of all q such that [i ? (C2g) : F] = 2 (which in particular requires 
that F contains the totally real subfield Q(C2<?) + of Q(C 2g )). Therefore, it 
suffices to compute the usual Dirichlet series or Euler product expansion 
for Cf(2) with the required precision; see also Dokchitser [9]. 

We now relate isometric circles to the arithmetic of B. Let p G ft have 
T p = {1}. A short calculation with the maps defined in (jl.ip shows that if 

^ G SL 2 (R), then = <$>g<$r x G SU(1, 1) has radius 

,/w 2Im(p) 
rad(J(g 9> )) = , 

/p(p) 
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where f g (t) = ct 2 + (d — a)t — b, a polynomial whose roots are the fixed 
points of g in C. We will abbreviate rad(g) = rad(7 (g^)) ■ The map 

. s invrad : M 2 (R) -► R 

1 j 5 ~|/*(p)| 2 + 2j/ 2 det(<,) 

yields a quadratic form on M2(M): explicitly, if p = x + yi, we have 
invrad ^ = (xa + b — (x 2 — ?/ 2 )c — xd) 2 + y 2 (a — 2xc — d) 2 + 2y 2 (ad - 

= y 2 (a — xc) 2 + (xa + 6 — x 2 c — x<f) 2 + y 4 c 2 + y 2 (xc + d) 2 , 

and hence the form invrad is positive definite and via too induces a positive 
definite form invrad : B — > R. For g & B, we note that dettoo(<7) = 
v(nrd(g)), where t> is the unique split real place of 5. 

Suppose that p = i. Then we have simply invrad ^ ^ = a 2 + b 2 + c 2 + 

d 2 . Let B = ( j . Identify F with its image F R under the unique 

split real place of B; without loss of generality, we may assume that h > 0. 
We may therefore embed : B ^ M 2 (R) by letting 

where sgn denotes the sign. Therefore if g = x + ya + z(3 + wa(3 G 5, then 
we see directly that 

invrad(g) = x 2 + hy 2 + \k\z 2 + h\k\w 2 . 

For the ramified real places v of F, corresponding to B <^-> i? (gip R = H, 
the reduced norm form nrd„ : -B — > R by g -u(nrd(g)) is positive definite. 
Putting these together, we find that the absolute reduced norm 

N : B -> R 

9 ^ 2 y 2 J2veS,v\oo nTd v(g) + fnvrad(#) = \f g {p)\ 2 + 2y 2 Tr F/Q nrd(p) 

is positive definite and gives the structure of a lattice of rank An. 

The elements g £ O with small absolute reduced norm N are those 
such that \ f g (p)\ and Trp/Q nrd(<?) are both small — in particular, this will 
include the elements of 0\ with small invrad (with respect to p G f)), which 
correspond to elements g £ T whose isometric circle in D (centered at p) 
has large radius. Since the Dirichlet domain D(p) has only finitely many 
sides, those g £ T with rad(g) sufficiently small radius cannot contribute 
to the boundary of D(p). 

Hence, one simple idea to construct D(p) would be to enumerate all 
elements of 0\ by increasing absolute reduced norm N until the exterior 
domain of these elements has area equal to /j,(T\$)). This method shows 
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that D{p) is indeed computable, and may have been known to Klein; it 
is mentioned by Katok [T7] when F = Q and sees further explication by 
Johansson [15]. Using the above framework, we can immediately improve 
upon this method by enumerating such elements efficiently using lattice 
reduction, as follows. 

Algorithm 3.2. Let OcBbea quaternion order. This algorithm returns 
a Dirichlet domain for T(O). 

1. Compute A = fj,(T(0)\Sj) by Remark EH 

2. Embed O w IR 4n as a lattice using the absolute reduced norm form 
N, and choose C G M>o. 

3. Using the Fincke-Pohst algorithm [10J, compute the set 

G{C) = {iooig/u) :±geO, N(g) < C, md(g) = u 2 G Zf} c T. 

4. From Algorithm 12.51 compute E = ext(G(C)). If /i(-E) = ^4 < oo, 
then return E; otherwise, increase C and return to Step 2. 

Remark 3.3. In choosing C, we note that 

{g G O; : rad( 5 ) > R} = L G O : iV( 5 ) < 2y 2 (» + ^)Wj; 

in practice, we would like to take C large enough so that G(C) / but 
not too large. It is not immediately clear how to choose C (and a strategy 
for its incrementation) optimally in general, unless one knows something 
about the radii of the sides of the Dirichlet domain. 

Our final algorithm (Algorithm 14.8ft significantly improves on Algorithm 
13.21 by the use of a reduction algorithm, which we introduce in the next 
section. 

4. Reduction algorithm 

In this section, we introduce the reduction algorithm (Algorithm 14. 3ft 
which forms the heart of the paper. This algorithm will allow us to find a 
normalized basis for the group T (Algorithm 14.7ft , yielding a fundamental 
domain. 

Throughout this section, let G = {gi, . . . ,gt} C T \ {1} be an (ordered) 
finite subset of a Fuchsian group T, and denote by (G) the group generated 
by G. For any z G D, we have a map 

p ■ r — M> 

7 i— ► p(-y; z) = d^z, 0) 
where d denotes hyperbolic distance. We abbreviate p(7;0) = p(7). 

Definition 4.1. Let z G D. An element 7 G T is (G, z)-reduced if for all 
g G G, we have ^(7; 2) < p(gj; z), and 7 is G -reduced if it is (G, 0)-reduced. 
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Remark 4.2. By Proposition 11.31 we note that 7 is (G, z)-reduced if and 
only if 72; G ext(G). 

We arrive at the following straightfoward algorithm to perform (G, z)- 
reduction. 

Algorithm 4.3. Let 7 G T and let z G 2). This algorithm returns elements 
7 G T and 5 G (G) such that 7 is (G, z)-reduced and 7 = £7. 

1. Initialize 7 := 7 and 5 := 1. 

2. If ^(7; z) < p{g 7 y; z) for all g G G, return 7, £. Otherwise, let g G G 
be the first element in G such that 

P(gr,z) = mm p(9iT,z). 
1 

Let 7 := <7i7 and 5 := and return to Step 2. 

We denote the output of the above algorithm 7 = rede (7; z) and abbre- 
viate redc(7;0) = redc(7)- 

Proof of correctness. The output of the algorithm 7 is by definition G- 
reduced. The algorithm terminates because if 7 1 ,7 2 , . . . are the elements 
that arise in the iteration of Step 2, then p(^)x, z) > p(7 2 ; z) > . . . ; however, 
the action of T is discrete, so among the points {7j(z)}i, only finitely many 
are distinct. □ 

A priori, Step 2 in Algorithm S3] depends on the ordering of the set G and 
therefore the output 7 will depend on this ordering. This is analogous to 
the situation of the reduction theory of polynomials, as follows. Let A; be a 
field, let R = k[x±, . . . , x n ] be the polynomial ring over k in n variables with 
a choice of term order, and let G = gi, . . . ,gt G R be not all zero. Applying 
the generalized division algorithm, one can reduce a polynomial / G R with 
respect to G, and the result is unique (i.e., independent of the ordering 
of the gi) for all / if G is a Grobner basis of the ideal / = (g±, ... ,gt). 
Moreover, if G is a Grobner basis, then / G / if and only if the remainder 
on division of / by G is zero. (See e.g. Cox-Little-O'Shea [SJ Chapter 2].) 
We can prove analogous statements, replacing the ring R by the group T, 
as follows. 

Proposition 4.4. Suppose that ext(G) is a fundamental domain for (G). 
Then for almost all z G D, redci'y^z) as an element of T is independent 
of the ordering of G for all 7 G (G). Moreover, for all 7 G T, we have 
redc(7) = 1 an d or ^V 1 £ (G). 

Here, "almost all" means for all z outside of a set of measure zero: it 
suffices to take z in the T-orbit of the interior of ext(G). 

Proof. Suppose that ext(G) is a fundamental domain for (G). Let z be in 
the T-orbit of zq G int(ext(G)), let 7 G (G), and let 7 = rede (7; z). Then 
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by Remark 14,21 we have jz E ext(G), and since ext(G) is a fundamental 
domain and Tz = Tzo with zq E int(ext(G)), we must have 72: = zq; in 
particular, 7 is unique and independent of the ordering of G. The second 
statement follows similarly: we have that 0,7(0) E int(ext(G)), so if 7 7^ 1 
then 7 0(G). □ 

Inspired by the preceding proposition, we make the following definition. 

Definition 4.5. A set G is a basis for T if ext(G) is a fundamental domain 
for (G) = r. If G is a basis that forms a normalized boundary for T, then 
we say that G is a normalized basis. 

Remark 4.6. It follows from Proposition 14.41 that if one can compute a 
normalized basis G for T, then one also has a solution to the word problem: 
given any element 7 E T, we compute 7 = rede 7, which by Proposition 
14.41 must satisfy 7 = 1, so we have explicitly written 7 as a word from G. 

We construct a normalized basis for as follows. 

Algorithm 4.7. Let G C T. This algorithm returns a normalized basis 
for (G). 

1. Let G := {g 1 , . . . , g t , g^ 1 , . . . , g^ 1 }. 

2. Compute the normalized boundary U of ext(G) by Algorithm 12.51 

3. Let G' := U. For each g E G, compute 5 = redc\{ ff } (g) using Algo- 
rithm H3J If 5 + 1, set G' := G' U {5}. 

4. Compute the normalized boundary J7' of ext(G'). If U' = U, set 
G := G' and proceed to Step 5; otherwise set U := U' and return to 
Step 3. 

5. If all vertices of E = ext(U) are paired, return U. Otherwise, for 
each g E G with a vertex t> E 1(g) which is not paired, compute 
~g := ied.G(g',v), where if v is a vertex at infinity we replace v by 
a nearby point in I(g~ l ) \ E C Q. Add the reductions g for each 
nonpaired vertex v to G and return to Step 2. 

Proof of correctness. First, note that if v be a vertex of E = ext(G), then 
by Corollary 11.41 v is a paired vertex if and only if for every side s C 1(g) 
containing v, we have that gv E I(g~ l ) is a vertex of 

Next, we prove that if the algorithm terminates it does so correctly. We 
construct a side pairing as in §1. A side s of E pairs up with gs C I(g~ l ) if 
and only if its vertices are paired, necessarily with the vertices of I(g~ l ) by 
Corollary 11.41 Therefore if we terminate in Step 5, we have in fact paired 
all sides of ext(t7) and by Theorem II . H ext(t7) is a Dirichlet domain and 
U is a basis. 

Otherwise, by Step 5 we have v E s such that gv ext(G). We now 
compute ~g = rede* (g; v), and refer to Proposition 11.31 Since v E 1(g), we 
have d(v,0) = d(gv,0), and since gv ext(G), we have d(gv,0) > d(~gv,Q). 
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Putting these together, we find that v G int(J(g)) and hence ext(GU{g}) C 
ext(G). 

Consider now the limit of the sets Goo = hmG and = limU as we 
let the algorithm run forever. Accordingly, every vertex v of ext(L r 00 ) must 
be paired, otherwise it would be caught in some step of the algorithm. 
Therefore by the above, Uoo is a basis for (Goo). But at each step of the 
algorithm, the group (G) remains the same, even as G changes: indeed, 
in Step 3, if g = 1 then g G (G \ {<?}). Therefore (Goo) = (G), and since 
(G) is finitely generated we know that U is finite, and hence the algorithm 
terminates after finitely many steps. □ 

We now extend this in the natural way to an arithmetic Fuchsian group 

no). 

Algorithm 4.8. Let O be a quaternion order. This algorithm returns a 
basis G for V = T(0). 

1. Choose G G R>o, initialize G := 0, and compute A = ii(T\$)). 

2. Using Steps 1-2 in Algorithm 13.21 compute the set G(G) C V. 

3. Call Algorithm |47l with input G U G(G) and let G be the output. If 
/i(ext(G)) = A < oo, then return G; otherwise, increase G and return 
to Step 2. 

A fundamental domain for an arithmetic Fuchsian group T C T(O) can 
easily be computed from this by first running Algorithm 14 . 81 and then com- 
puting a coset decomposition of T in T(0); and for that reason, one may 
even restrict consideration to the case where O is maximal. 

Remark 4.9. In practice, in some cases we can improve Step 5 of Algorithm 
14.71 for arithmetic Fuchsian groups as follows. For each nonpaired vertex 
v, we can consider those elements with small absolute reduced norm N 
relative to p G D taken to be a point along the geodesic between and 
v: indeed, by continuity if g G 0\ has v G int(/(g)), then rad(g) increases 
as the center p moves towards v and thus N(g) decreases, so using lattice 
reduction we are likely to find a small such g. 

5. Proof of the main theorem 

We are now ready to prove the main theorem of this paper. 

Theorem 5.1. There exists an algorithm which, given a finitely generated 
Fuchsian group T and a point p G with T p = {1}, returns the Dirichlet 
domain D{p), a side pairing for D{p), and a finite presentation for V with 
a minimal set of generators. 

To prove the theorem, we need to show how the output of Algorithm 14.71 
yields a finite presentation for V with a minimal set of generators. Indeed, 
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Algorithm 14,71 terminates only if it has computed a side pairing P (which 
we may assume meets the convention in §1) for the Dirichlet domain D. 
Such a side pairing P gives a set G of generators for V by Proposition 11.11 

We now consider the induced relation on the set of vertices. A cycle 
of D is a sequence vi,...,v n = V\ which is the (ordered) intersection of 
the T-orbit of v = v\ with D. To each cycle, we associate the word g = 
9n9n-i ■ ■ ■ 929i where gi(v{) = Vi + \ and the indices are taken modulo n. We 
say that a cycle is a pairing cycle if <?j G G for all i, and without further 
mention we shall assume from now on that a cycle is a pairing cycle. 

A cycle is minimal if Vi ^ Vj for all % ^ j. Every vertex v of D is contained 
in a unique minimal cycle (up to reversion and cyclic permutation). Indeed, 
by the uniqueness of the side pairing, a vertex v G 1(g) D I(g') either 
has v = gv = g'v, in which case v has nontrivial stabilizer and one has 
the singleton cycle v, or v has trivial stabilizer and is paired with the 
distinct elements gv G I(g~ l ) and g'v G /(g' -1 ), each of which also has 
trivial stabilizer, and then continuing in this way one constructs a (unique 
minimal) cycle. This analysis gives rise to the following algorithm. 

Algorithm 5.2. Let P be a side pairing for a Dirichlet domain D for T. 
This algorithm returns a set of minimal cycles for D. 

1. Initialize V to be the set of vertices of D and M := 0. 

2. If V = 0, terminate. Otherwise, choose v G V with w = /(g) D /(</) 
for g,g' G G(P). If gv = t>, add the cycle v to M and remove u from 
V, and return to Step 2. Otherwise, let i := 1 and v\ := v. 

3. Let fj+i := gvi G I(g~ l ) H /(</)• ^ = u l> a dd the cycle . . . , 
to M, remove these elements from V, and return to Step 2; otherwise, 
increment z := i + 1, let 5 := (/ and return to Step 3. 

The relations associated to minimal cycles have the following important 
property. 

Lemma 5.3. Let g G G be a side-pairing element. Then g appears at most 
once in any word associated to a minimal cycle. Moreover, g and its inverse 
appears in exactly two such words. 

Proof. By definition, a side-pairing element g pairs a unique set of sides: 
in particular, g pairs the vertices of one side s with the vertices of another. 
Suppose that g occurs twice in a word associated to a minimal cycle. Then 
by minimality, the vertices of s are in the same T orbit. But this implies 
that g maps 1(g) to itself, so g has order 2 and therefore one of the vertices 
of s is fixed by g, a contradiction. 

In a similar way, we see that g and its inverse can appear in at most two 
words since each vertex belongs to exactly one minimal cycle. □ 

We have the following characterization of the minimal cycles. 
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Proposition 5.4. Beardon [21 Theorem 9.4.5] For all p £ $j outside of a 
set of area zero, the following statements hold: 

(i) Every elliptic cycle has length 1; 

(ii) Every accidental cycle has length 3; and 

(iii) Every parabolic cycle has length 1. 

Remark 5.5. The exceptional set of p is contained in the union 

E 2 = \J {z : R{z) G M} 

over all triples f,g,h S T such that 

\ = {z~gz){fz - hz) 
{) ( z -fz)(gz-hz) 
is not constant. It is easy to see that the set E 2 has area zero. 

For the purposes of computing a minimal set of generators and relations, 
we may and do assume that p does not lie in the exceptional set; indeed, 
a sufficiently general choice of p will suffice, and so in practice the condi- 
tions of Proposition 15.41 always hold. In particular, every elliptic cycle is 
represented by a minimal cycle (whose fixed point is a vertex of D). 

Now, to each cycle, associated to the word g, we further associate a 
relation in T as follows. By definition, we have g € T v , and therefore we 
have one of three possibilities. If #T V = 1, then we have the relation g = 1; 
we call g an accidental cycle. If 1 < #T V < oo, then we associate the 
relation g k = 1 where k is the order of g, and we call g an elliptic cycle. 
Otherwise, if #T V = oo, then we associate the empty relation, a parabolic 
cycle. We note that the latter occurs if and only if g has infinite order if 
and only if trd(g) = ±2, so the relation g is computable. 

We now appeal to the structure theory for Fuchsian groups with cofi- 
nite area [TBJ §4-3]. Suppose that V has exactly t elliptic cycles of orders 
mi,..., mt £ <^>2 and s parabolic cycles, and that X = T\f) has genus 
g. We say then that T has signature (g; mi, . . . , mt; s). Moreover, T is 
generated by elements 

(5.1) cti,... ,a g ,Pi,.. . ,/3 9 ,7i,. . . ,7t,7t+l> • • -,Jt+s 
subject to the relations 

(5.2) 7i mi = • • • = 7t mt = K A] • • • [<* g , P 9 ]ll ■ ■ ■ lt+s = 1, 

where [a, f3] = afia -1 is the commutator. (One obtains a minimal set of 
generators from this presentation by eliminating jt+s whenever t + s > 0.) 

From the set of generators coming from the side-pairing elements and the 
set of relations coming from the minimal cycles, we can build a minimal 
set of generators and relations by "back substitution". First, we prove a 
lemma. 
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Lemma 5.6. Suppose V = Ti *T2 is a free product, and that 73 G Ti or T2 
for i = 1, . . . , s + t. Then either Ti or T2 is isomorphic to the free product 
of cyclic groups. 

Proof. Let <f> : T — > Ti * r 2 be an isomorphism. Passing to the quotient by 
the 74, for i = 1, . . . , t + s, we may assume that s = t = 0. But then the 
homology groups Hi(T, Z) (coming from group homology) coincide with the 
homology groups Hj(Y, Z) (coming from topology) where 1" is the orientable 
surface of genus g [S §11.4]; in particular, we have Hq(T,Z) = Z. By the 
Mayer- Vietoris sequence jSJ Corollary II. 7.7], we have 

z * h (t, z) * # (ri * r 2 , z) flb(ri, z) e H (r 2 , z) 

so say i?o(r2, Z) = 0; but this immediately implies T2 is trivial as well, and 
the result now follows. □ 

Algorithm 5.7. Let P be a side pairing for D and let M be a set of 
minimal cycles for D. This algorithm returns a minimal set of generators 
and relations for T. 

1. Let -ff C G(P) be such that g G G implies either g = g -1 or g -1 G. 

2. Let -R be the set of elliptic cycles in M and let A be the set of 
accidental cycles. Initialize r to be an element of A and remove r 
from A. 

3. If A = 0, add r to R and return the generators H and the relations i?. 
Otherwise, choose an element g £ A such that g and r have an element 
gi £ H in common; then solve for gi, substitute this expression in for 
gi in the relation r, and remove gi from H. Return to Step 3. 

Proof of correctness. If in Step 3 there is always an element g G A such 
that g and r have an element in common, then the algorithm terminates 
correctly: in the notation of (|5.1H5.2|) . there are exactly t + l relations, and 
hence the set of generators must also be minimal. 

So suppose otherwise. Let H\ be the set of g G H such that g or g~ l 
occurs in the relation r and let H2 = H \ H\. Let Ti,r2 be the groups 
generated by H\,H2- Then by assumption, V is the free product of Ti and 
T2- By Lemma 15.61 since the relation in T\ is nontrivial, it follows that I^ 
is the free product of finite cyclic groups, and hence cannot contain any 
accidental cycles, which is a contradiction. □ 

The minimal presentation resulting from Algorithm 15.71 is not necessar- 
ily of the form (|5.1[) — f)5.2[) ; we refer to the methods of Imbert [13] for an 
alternative approach using fat graphs which computes such a canonical 
presentation. 

This completes the proof of the theorem and the accompanying corollar- 
ies in the introduction. 
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Remark 5.8. If in the first corollary, one wants the structure of O* , we use 
the exact sequence 

1 _> if o* -> o* ^ z* Ft+ /zf -» 1 

where Z*^ + = {u £ : u(u) > for all ramified places v \ oo}. From 
the solution to the word problem, it then suffices to find elements 7 S O* 
such that nrd(7) = u generates the finite group r L* F+ /'L*f : and these can 
be found using the methods of §3. 

6. Examples 

We have implemented a variant of the above algorithm in the computer 
system Magma [3]. In this section, we provide two examples of the output 
of this algorithm. 



1 

Figure Ol: A Dirichlet domain for the arithmetic Fuchsian group Fg(13) 

/3, — 1\ 

First, we consider the quaternion algebra B = I — ^— — I of discriminant 



6. A maximal order O is given by 

o = z®Zaezp® I } + a + P + a P _ 

We consider the Eichler order contained in O of level 13, given by 

0(13) = Z0Z 3 ~ 5a ~ 2 5/3 + 3a/3 @ Z(2-2a-/3 + a/3) 
13- 13a - 13/3 + 13a/? 
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We denote T(0) = rg(13). We embed B ^ M 2 (R) by the embedding 
(|3.4p . and take p = 9i/10 G f}. By (|3,ip . we compute that the Fuchsian 
group I(j(13) has coarea 14/3. 

Step 2 in Algorithm 14, 81 finds the units (1— a— 3/3+a/3)/2, a— 2/3, . . . , and 
following the algorithm, reduction and further enumeration automatically 
yields the fundamental domain as in Figure [BJl. (The methods in Magma 
for producing the postscript graphic are due to Helena Verrill |21j.) 

This domain already exhibits significant complexity: it has 38 sides and 
hence 19 side-pairing elements, which yields a set of 10 minimal generators 
71, . . . , 710 for r(j(13), namely 

12 - 7a + 4/3 + 2a(3, (1 - a - 33/? - 19a/3)/2, 2a + 16/? + 9a/3, 
(37 - 19a + 9/3 + lla/3)/2, 2a + 4/3 + a/3, (1 - a - 3/3 + a/3)/2, 
a -2/3, (l + 7a- 15/3 - 5a/3)/2, (1 + 7a - 45/3 - 25a/3)/2, a - 14/3 - 8a/3, 

subject to the relations 

73 = 7 5 2 = 7? = 7i 2 o = 7 2 3 = 7e 3 = 7 8 3 = 7 9 3 = 1 
7f 1 74757fr 1 7i72~ 1 737^ 1 777^ 1 7 9 _1 7ro 1 = 1. 

We deduce that Tq(13) has signature (1; 2, 2, 2, 2, 3, 3, 3, 3; 0), a fact which 
can be independently verified by well-known formulae pQ. 

Second, we consider the totally real number field F generated by a root 
t of the polynomial x 7 — x 6 — 6x 5 + 4x 4 + 10x 3 — 4x 2 — 4x + 1; it is the 
minimal septic totally real field, having discriminant dp = 20134393 = 
71 • 283583. We consider the quaternion algebra B which is ramified at 6 



of the 7 real places of F and no finite place: explicitly, B = I — — I where 



h = -t 6 + 6t 4 + 1 3 - 9t 2 - 3t + 1 and k = -t 2 + 2t - 1, and in fact h, k G Z* F . 
We compute a maximal order O of B. Letting T = r(O), we see that T 
has coarea 5/2. The output of Algorithm 14.81 in this case is given in Figure 
EJ2; we find that V has signature (0; 2, 2, 2, 2, 2, 3, 3, 3; 0). 

We conclude by noting that it would be interesting to extend the methods 
in this paper to other arithmetic groups; this would allow the computation 
of unit groups for a wider range of quaternion algebras over number fields 
and would have further consequences for the algorithmic theory of Shimura 
varieties. 
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